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Inconel 718 is a nickle-based superalloy that possesses 
several outstanding elevated- temperature mechanical properties. It has been 
widely used for the manufacture of critical compressor and turbine engine 
components requiring relatively long service lives. A large number of 
tests, including tensile, creep, fatigue, and creep-fatigue have been 
performed to characterize the mechanical properties of Inconel 718 at 
1200°F, the operating temperature for turbine blades. In addition, a few 
attempts have been made to model the behavior of Inconel 718 at 1200 F using 
viscoplastic theories. 

The Chaboche theory of viscoplasticity can model a wide variety of 
mechanical behavior, including monotonic, sustained, and cyclic responses of 
homogeneous, initially- isotropic , strain-hardening (or softening) materials. 
It has been successfully used to model the viscoplastic behavior of several 
structural materials of practical importance. This paper shows how the 
Chaboche theory can be used to model the viscoplastic behavior of Inconel 
718 at 1200°F. First, an algorithm has been developed to systematically 
determine the material parameters of the Chaboche theory from uniaxial 
tensile, creep, and cyclic data. The algorithm, however, is general and can 
be used in conjunction with similar high- temperature materials. A 
sensitivity study was then performed and an 'optimal' set of Chaboche 's 
parameters was obtained. This study has also indicated the role of each 
parameter in modeling the response to different loading conditions. Based 
on the 'optimal' set of material parameters, uniaxial tensile, creep, and 
cyclic behavior has been predicted. The results were compared to available 
experimental data and the agreement was found to be good. Moreover, 
predicted behavior of Inconel 718 at 1200°F under a variety of additional 
loading conditions has been examined and relevant conclusions were drawn. 

I . INTRODUCTION 


For high temperature applications, such as those encountered in the 
nuclear and aerospace industries, severe demands are generally placed on 
candidate structural materials. These demands have led to the development 
of a new class of structural materials called superalloys having 
considerably improved elevated- temperature mechanical properties. One such 
material is Inconel 718, which possesses several excellent properties such 
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as creep -rupture strength, high-cycle fatigue strength, oxidation 

resistance, and long time stability [1-3]. These properties are, for the 
most part, due to the susceptibility of this material to precipitation 
hardening and strengthening caused by coherent gamma prime precipitates. 
These properties can be obtained by two different solution annealing and 
aging procedures. A comparative study of the effects of both types of heat 
treatment procedures on the tensile and fatigue behavior of Inconel 718 from 
room temperature to 1200°F has been presented in Ref. [2], It was concluded 
that Inconel 718 subjected to the so called 'standard heat treatment' 
exhibits higher resistance to stress rupture and fatigue. Therefore, 
Inconel 718 in this condition is widely used in current production of gas 
turbine engine components requiring relatively long service lives and 
operating normally at high temperatures [1,4]. 

A large number of uniaxial monotonic, sustained, and cyclic tests have 
been performed to obtain the mechanical properties of Inconel 718 at 1200°F 
[1-9]. However, only a few of these tests address the strain-rate 

sensitivity, with the result that no significant strain rate sensitivity 

could be observed for e > 5x10 ^ sec \ At lower strain rates, tensile and 

cyclic data were too limited to establish the rate sensitivity of this 

material, although existing creep data do suggest the existence of some rate 
sensitivity. 

Domas, et al [4] have shown that the tensile and stabilized cyclic 
stress -inelastic strain curves for Inconel 718 at 1200°F can be described by 
the Ramberg-Osgood form 

* n* 

a - K (€") n , (1) 

If 

where a and c are the tensile stress and inelastic strain (or the stress at 
the tensile tip of the stabilized hysteresis loop and the corresponding 
inelastic strain amplitude), respectively, and K and n are material 

parameters that depend on the heat treatment and strain rate. Values of K 

^ o 

and n for Inconel 718 (standard heat treatment) at 1200 F and for different 

tensile and cyclic strain rates are listed in Table 1. The first five sets 

of values were determined by Domas, et al [4], whereas the remaining values 

were determined within a larger study [10], part of which is the present 
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work. Available experimental data on the cyclic behavior of Inconel 718 at 
1200°F have also indicated that this material cyclically softens for 
approximately 10-20% of its fatigue life, followed by a period of near- 
stabilization when damage mechanisms are activated. For convenience the 
values of a and e used to evaluate the cyclic material parameters in Table 
1 were those at N^/2, where is the fatigue life. 

Inconel 718 at 1200°F displays classical creep response (primary, 
secondary, and tertiary creep stages). The time to the onset of tertiary 
creep is normally about 0.5 t r> where t is the rupture life [11]. For 
Inconel 718, however, the region of increasing creep rates following steady 
state creep commences at about 0.3 t [3]. This suggests that the 
initiation of increasing creep rates may be due to softening and is not 
necessarily a manifestation of the instability associated with tertiary 
creep [3,6]. Thus, from a unified viscoplastic theory view -point, some 
provision for softening (or increasing creep rates) must be incorporated 
into the theory. Booker [3,6] has developed a creep model in the form 

€ - EXP [1.75(t - 1)] • t°' 2 , 0 < t < 1, (2) 

n n n n 

which describes the salient features of the creep behavior of Inconel 718 at 

1200°F. In Eq. (2), e - e/e is the normalized creep strain, t - t/t 

^ n tr n tr 

is the normalized time; provided that € and t are the strain and time to 
tertiary creep given empirically by 


- 0.392 t 1 * 04 , 

r 

(3) 

- 0.2 + 1.16 t" 0,14 • 

(4) 


The time to failure, t , is given by 


log t - 163.92 + 23924/T - 226.92 log <? h 

+ 94.196 (log a) 2 - 13.215 (log a) 3 , 


(5) 
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where T is the temperature in °K and is the hold stress in MPa. In Eqs . 
(2) thru (5), time is given in hours and strain in percent. 

A few attempts have been made to model the behavior of Inconel 718 at 
1200°F using unified viscoplastic theories, e.g., [7,8]. In this paper, the 
Chaboche yield-based theory of viscoplasticity [12,13] is used for the same 
purpose. This theory can model a wide variety of rate -dependent inelastic 
mechanical behavior, including monotonic, sustained, and cyclic responses. 
It also models isotropic and kinematic strain hardening, including cyclic 
hardening (or softening) and the multidimensional Bauschinger effect. This 
theory has been employed to model the mechanical response of Inconel 100 
[12,14], type 316 stainless steel [15,16], and Ti-6Al-4V alloy [17,18]. 
However, the theory does not account for the rate dependence of initial 
yield, which limits its use to materials that do not exhibit such behavior 
within the loading rates of interest. In a previous paper [19], the authors 
have shown how rate dependency of initial yielding can be incorporated into 
the Chaboche theory, without changing its general structure. The principal 
objective of the current paper is to show how the extended theory can be 
used to model the mechanical behavior of Inconel 718 at 1200°F. A major 
part of this effort has been devoted to developing a systematic algorithm 
for the determination of the Chaboche material parameters from available 
uniaxial test data. The need for a systematic procedure for determining the 
material parameters of a viscoplastic theory from the results of standard 
mechanical tests has been addressed [17,21]. 

A sensitivity study was then undertaken to explore the effects of 
introducing small changes in the Chaboche material parameters on different 
predictions and to obtain an 'optimal' set of the material parameters that 
globally improves the predicted behavior. The 'optimal' set was then 
incorporated into the Chaboche set of differential equations; and uniaxial 
tensile, creep, and cyclic responses were predicted and compared to 
available experimental data. In addition, strain- rate and strain- rate - 
history effects, creep, stress relaxation, and the load- unload- reload 
effects were predicted and examined. 
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II. CONSTITUTIVE EQUATIONS 


In the present paper, the one -dimensional form of the Chaboche theory 
is presented and used; the three-dimensional form of the theory can be found 
elsewhere [12,13]. A basic assumption of the theory is the decomposition of 
the rate of deformation into rate- independent elastic and rate -dependent 
inelastic components. For infinitesimal deformations, the strain rate, e, 
approximates the rate of deformation, or 

l - c' + l" • (6) 

The linear elastic strain rate, l 9 , is obtained from the time derivative of 
Hooke's law, that is, 

e ' — a/E (7) 


where a is the stress and E is Young's modulus. The nonlinear inelastic 

• " 

strain rate, e , is based on the normality hypothesis and has the form 


r 11 . 


KJ Sg(<7 - Y) , F > 0, 
0, F < 0, 


( 8 ) 


where K and n are material parameters. F is the von Mises yield function 
given by 

F(<r,Y,p) - jo - Y| - R(p) , (9) 


where Y is the kinematic hardening variable, p is the cumulative inelastic 
strain defined in terms of the inelastic strain rate by 



( 10 ) 


and R is the isotropic hardening variable associated with p. Note that in 
Chaboche 's theory, as in most of the viscoplastic theories, two variables, Y 
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and p, are employed to model hardening. Note also that F represents the 
overstress, whose positive values signify the initiation of inelastic 
deformations [c.f. Eq. (8)]. 

The evolution of the kinematic hardening variable is given by 

Y - c(ae " - Yp) - 7 |Y| m Sg(Y) , (11) 

where c,a, 7 , and m are material parameters. The first term on the right 

hand side of Eq. (11) is the nonlinear Armstrong -Frederick modification of 

Prager's linear kinematic hardening rule, which is appropriate for cyclic 

loading. The second term models softening effects such as reduced hardening 

rate, secondary creep, and stress relaxation typically associated with 

elevated temperatures. Equation (11) provides an effective tool by which 

anisotropic hardening (Bauschinger effect) and creep/relaxation effects can 

be appropriately modeled and represents one of the principal advantages of 

the Chaboche theory. In its modified form [19], the isotropic hardening 

variable, R, depends on the cumulative inelastic strain, p, and the total 

strain rate, e, rather than on p alone, as postulated in the original 

theory. The dependence of R on c is assumed to be limited to its initial 

value, R , that is, 
o 

R(p.O - q + [R q (0 - q] e" bp , (12) 

where b and q are material parameters and R q is the rate -dependent initial 
yield strength. 

III. DETERMINATION OF THE CHABOCHE MATERIAL PARAMETERS 

In its modified form, the Chaboche theory incorporates eight 
viscoplastic material parameters and the rate -dependent initial yield 
strength [c.f. Eqs. (8), (11), and (12)], which are generally temperature- 
dependent. It is now shown how these parameters as well as the function 
R Q (e) were determined from available experimental data of Inconel 718 at 
1200°F. Some approximations and assumptions had to be made; these will be 
discussed as they are introduced. 
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Determination of a and c 

The Chaboche flow law given by Eqs . (8) and (9), upon inversion, takes 
the form 


o - Y + R + 


K(J") 1/n • 


(13) 


At relatively high strain rates, the recovery term in Eq. (11) can be 
neglected, thereby allowing it to be integrated as 


where 


-sc(€ m - e") 

Y - sa + (Y - sa) e ° 


s - Sg(^") 



> 0 , 

e" < 0, 


(14) 


(15) 


and Y and € M are the initial values of Y and c". 
o o 

For the stabilized cycle of a strain-controlled cyclic test, R assumes 

its saturation value q [c.f Eq. (12)], whereas Y alternates between two 

extreme values, Y and Y . - - Y . It can be shown that [10,12] 

max min max J 

- a U - e' CAe ")/(l + e' cAe ”) - a tanh(ce") , (16) 

max a 

where - Ae"/2 is the stabilized inelastic strain amplitude. Hence, if a 
series of stain-controlled cyclic tests are performed until stabilization 
has occurred, a cyclic stress -strain curve can be established through the 
tensile peaks of these stabilized hysteresis loops. Such a curve, within 
the context of Chaboche' s theory, can be represented by use of Eqs. (13) and 
(16) as 

a - a tanh(cc w ) + q + K(«")^ n , (17) 

a 

which, when differentiated with respect to 6^, yields 

do/de - a c sech^(c€ w ) • (18) 

a a 
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An alternative form for Eq. (18) can also be obtained by 
differentiating Eq. (1), i.e., 

* 

« * * n - 1 

da/de - n K («") • (19) 

a a 

Since values of K and n are available (c.f. Table 1), Eq . (19) can be used 

to generate values of f. - do/de M for (c 11 )., i— l,2,...,n. These values can 

be used to evaluate a and c in Eq. (18). This is, in fact, a nonlinear 
regression analysis problem. Draper and Smith [22] have shown how such a 
problem can be solved based on the least squares technique. A computer 
program based on their analysis has been developed and used to find the 

values of a and c listed in Table 1. The first and seventh sets are 

considerably different from the remaining values. Therefore, these two sets 
were excluded and mean values of a and c were calculated for the remaining 
sets as 

a - 31.40 KSI , c - 349.5 • (20) 

Determination of K and n 

Since the material parameters K and n characterize the rate dependency 
of the material response, they were evaluated from tensile data at different 
strain rates. Two tensile stress-strain curves may be adequate for this 

purpose. Figure 1 illustrates two stress-strain curves at strain rates of 

• • • * 

and € where If it is assumed that these rates are 

sufficiently high and close together, then Eq. (14) holds for Y, and R can 

o 

be considered the same for both cases. Thus, 


Y - Y , Y. - Y. , R - R 
a l a 2 b l b 2 a l a 2 


vs 


( 21 ) 


where a^,a 2 ,b^, and are shown in Fig. 2. 

From Eqs. (13) and (21), the stress differential may be written as 


ho - o 
a a r 


- o - K |(e" 

a x L a 


)Vn . ( - ,1/n _ 


( 22 ) 
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which relates the difference in stress levels to the corresponding inelastic 
strain rates at specific inelastic strain values. It is necessary, 
therefore, to evaluate the inelastic strain rate at different points on both 
curves using the equation 

e" _ l . \< - e - a/E - e - J(^)/E - «(1 - E*/E) . (23) 

'fc , 

where E - da/de is the instantaneous slope of the stress -strain curve. 

Equation (1), along with the integrated form of Eq. (6), can be employed to 

obtain the relation 


_ 1 

da E 


+ 



(24) 


Note that the inverse of this equation is E . Therefore, Eqs. (23) and (24) 

can be used to calculate values of c" and € n in Eq. (22), and A a can be 

a l a 2 a 

obtained by applying Eq. (1) independently to curve 1 and curve 2 of Fig. 1. 
Then, K and n in Eq. (22) are evaluated using a nonlinear regression 
procedure similar to that used in evaluating a and c. The tensile data of 
Table 1 were used for this purpose, and the following values of n and K were 
obtained 

n - 5.12, K - 154.4 KSI n Jsec • (25) 


Determination of b and q 

Once a, c, n, and K have been obtained, b and q can be determined from 
the results of one strain-controlled cyclic test, such as that reported in 
Refs. [7,9] at l^l - 2.667x10 ^ sec ^ and c - ± 1%. At the tensile tip of a 
hysteresis loop, Eqs. (13) and (16) can be combined and solved for R as 

R - a - a tanh(ce") - K(J") 1/n . (26) 

a 


and p can be calculated from 


n 

p - e"/2 +22 Ae J , 

1 1-2 1 
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i 


( 27 ) 
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where n represents the n-th cycle and Ac? represents the inelastic strain 

range at the i-th cycle. Equations (26) and (27) were then used to 

calculate R and p values from the strain-controlled test discussed above, 

with the results given in Table 2. In these calculations the overstress, 

K(c" ) 1//n , was assumed to be constant, since the slopes of the loops at their 

peaks were relatively constant. The value of c" was calculated from Eq. 

-3-1 

(23) to be 2.465x10 sec , from which the overstress was determined to be 
47.73 KSI . Since the available results did not include records of all 
hysteresis loops, it was assumed that a missing loop had the same value of 
Ac” as the first preceding available loop. Table 2 shows that R had 
essentially saturated by the 32nd cycle. The corresponding value of R was, 
therefore, taken to be the saturation value, i.e., 


q - 18.63 KSI • (28) 

A value for b was then determined from the slope of the line 

ln(R-q) - ln(R Q -q) - bp, (29) 

to be 

b - 4.679 • (30) 


Determination of R 
o 

In a previous paper [ 19 ] , the authors have shown how the Chaboche 

theory can be modified to account for the rate -dependent initial yield 

response and permit prediction of creep behavior at low stress levels. In 

effect, the function R (e) in Eq. (12) had to be evaluated. Based on the 

o 

analogy between tensile and creep responses [23] and the creep data base of 
Ref. [8], a form for R (e) was selected as 


R 

o 


R 


oo 


r 

r 

r 

• 1 + ln-1 ] 

k 

L + EXP 


3 

2 

i-0 


^(lnO 



(31) 


where R^-q-18.63 KSI, c» o -25.56, Qj-4.551, q 2 -0.3360, and a 3 -0 . 8375x10' 2 , as 
shown in Ref. [19]. 
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Booker's creep model [c.f. Eqs. (2-5)] was used to evaluate the strain 
at the initiation of steady-state creep needed for these calculations. This 
was achieved by differentiating Eq. (2) twice with respect to time, equating 

to zero, solving for t^, and then substituting back into Eq. (2) to obtain 

the required strain values. 

Determination of m and -v 

It has been shown in Ref. [17] that hardening is mostly kinematic 

during primary creep, because of the relatively high value of the ratio of 

the hardening exponents c and b [c.f. Eqs. (12) and (13)], and because of 

the small amount of inelastic strains that are likely to have accumulated. 

Therefore, it can be assumed that R - R - constant at the attainment of 

o 

steady- state creep. During secondary creep, e ■ «" - constant, and the time 
derivative of Eq. (13) reduces to 

Y - c(a-Y) e gs - 7 Y° - 0, (32) 

e gg is the steady-state creep rate. Equation (32) can also be written 

ln[c(a-Y) e ] - m InY + lny, (33) 

s s 

which represents a straight line whose slope is m and intercept is In 7 . In 
this equation, Y can be calculated from Eq. (13) with R - R q . The data of 
Ref. [ 8 ] for small creep stresses, where the recovery effects are 
significant, were used to evaluate 7 and m as 

7 - 0.4x10 * 10 , m - 6.942 • (34) 

To conclude, a procedure has been developed to evaluate the material 
parameters of the modified Chaboche theory of viscoplasticity from uniaxial 
tensile, creep, and cyclic test data. Since these tests were not 
specifically designed for the purpose of determining the Chaboche 
parameters, the procedure employed herein to determine these parameters is 
not unique, but appropriate for the kind of data available. However, the 
procedure is general and can be employed to determine the Chaboche 


where 

as 
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parameters for any similar material. Values of the parameters are listed in 
Table 3. Based on these values, predictions of a tensile test at e - 
1.333x10 "^sec \ a creep test at 130 KSI, and a strain-controlled, fully- 
reversed cyclic test at | e | - 4x10 ^sec ^ and Ae - 2% are shown in Figs . 2 
to 4, respectively, along with available experimental data for comparison. 

IV. SENSITIVITY STUDY 

Although the predicted and actual responses are in reasonable 
agreement, it was anticipated that better agreement could be achieved by 
modifying the parameter values. Therefore, a parametric study was 
undertaken to explore the sensitivity of each of the parameters to the 
overall predicted behavior of Inconel 718 at 1200°F, and thus lead to the 
determination of an 'optimal' modeling of overall material response. To 
perform such a multi-parameter optimization, it is necessary to have a 
complete and accurate set of experimental results of specific base-line 
tests. These tests should be chosen in such a way that all aspects of 
material behavior to be modeled by the theory are represented in the test 
data. Such a data set is not available for Inconel 718 at 1200°F, partly 
because of the expense of such testing, but also because existing test data 
were not acquired for the purpose of evaluating Chaboche's model. 
Therefore, it has been necessary to work with available data sets that are 
not adequate for the present purpose. 

With this in mind, the data base of Ref. [8], which consists of tensile 
and creep test results only was employed to perform the parametric study. 
Then this data base was enlarged by including cyclic data from Ref. [4]. 
The first step of the parametric study consisted of varying the parameters, 
one at a time, by +7% of the value listed in Table 3, predicting the 
corresponding material responses, and comparing these predictions with the 
selected data base. Because of space limitations, only representative 
samples of the results of this step are shown in Tables 4 to 6 , along with 
corresponding experimental data. A thorough examination of the results of 
this step of the parameteric study has led to the following observations: 

(a) Some parameters have stronger effects on material behavior than 
others, in decreasing order, R q , n, k, and m have the strongest effects. 
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(b) Whereas a material parameter or a subset of the parameters has a 
strong effect on a certain facet of material behavior, it does not 
necessarily have the same effect on a different aspect of material behavior. 
For instance, R q had a strong effect on the creep behavior; however, it had 
less effect on the tensile and cyclic behavior. 

(c) Creep response seems to be the most sensitive to changes in 
parameter values. For instance, changing R q by ±7% results in approximately 
one order-of -magnitude change in the minimum creep rate value. 

(d) Increases in b, 7, and m increase the minimum creep rate and the 
inelastic strain range, and decrease the saturation stress, whereas the 
other remaining parameters have the opposite influence. 

Based on these observations, it was decided to extend the parametric 

study for the parameters that strongly affect the material behavior by 

reducing the amount of change from +7% to ±3.5% for n, K, and m and to 

+5.25, 3.5, and 1.75% for R , since it seemed to be the most influential. 

- o 

Representative results of these changes are also included in Tables 4 to 6. 

These results support the previous observations in that they confirm that 

the overall material behavior is most sensitive to variations of R . For 

o 

example, a change of +1.75% in R^ changes the minimum creep rate by about 50 
to 100%. Thus, the results of the current parametric study verify that the 
modification introduced to Chaboche's theory, i.e., considering R q to be a 
function of the strain rate, is of critical importance for improving the 
accuracy of this theory. 

Although the results of this parametric study have shown that the 
cyclic material response (first cycle Ac and Ac”) is not strongly affected 
by small changes in b and q, the influence of these parameters is best 
exhibited by studying the complete cyclic behavior and not just the first 
cycle. Unfortunately, the only cyclic test available for this material is 
the one previously used to determine values of b and q. This test was 
conducted at |c| - 2.667x10 sec , a value which is well above the limit 
of rate sensitivity for Inconel 718 at 1200°F - estimated to be about 4x10 
sec ^ . Thus the results of such a test should not be used to determine the 
material parameters. Problems also arise in attempting to use Chaboche's 
theory to predict a rate- insensitive response, since it is basically a rate- 
dependent plasticity theory. This problem was demonstrated by predicting 


49 


* 3 X 

the cyclic behavior at |e| - 2.667x10 sec and e - +1% for the first 

several cycles. The results showed that the stress was initially 
overstimated and the inelastic strain was underestimated. Also, the 
predicted cycles softened at a higher rate than actually observed, which 
suggests that b and q must be re-evaluated. 

To best accomplish such a task, the hysteresis loops from a cyclic test 
at | e | < 4x10 ^ sec ^ are needed. Such data, however, are not available*, 

and it was hypothesized by virtue of the material rate insensitivity for 
e > 4x10 sec that the results of a cyclic test at \e\ - 2.667x10* sec" 
are essentially the same as a test at |e| -4x10 ^ sec In effect, values 
of R have been recalculated and the new estimates of the material parameters 
b and q were found to be 


b - 3.75, q - 50 KSI • (35) 

It was also found that by rounding off the values of a, c, K, m, and n, 
i.e., 

a - 30 KSI, c -♦ 350, K -*■ 155 KSI n Jsec, m -*• 7, n - 5.1, (36) 


insignificant effects on the predicted behavior were observed. Therefore, 
these values will be considered in the remainder of this work to be the 
'optimal' values. 

The final step in the optimization procedure was to change the form for 
R q (0 such that the best possible agreement could be obtained between 
experimental and predicted behavior. First the values of R q that gave best 
prediction of overall behavior were determined and are shown in Ref. [19]. 
Then, an improved expression for R was established as 


R 


o 


R 


00 


r 

( 

r ^ 

« 1 + In i 

1 + EXP Ls n 

l 

V 

Li-0 


^ i (lne) 



(37) 


where R qq - q - 50 KSI, 0 q - 48.88, ^ - 18.30, ^ - 2.767, ^ - 0.2082, 

- 0.7813xl0* 2 , and /? 5 - 0.1174x10* 3 . 


Contacts were made with 0RNL to obtain the cyclic data from the tests 
reported in Refs. [2,5] but it was indicated that they were not available. 
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The revised material parameters of the extended Chaboche theory are 
listed In Table 3, and were employed to Improve the material response 
predictions In Figs. 2 to 4. 


IV. PREDICTIVE CAPABILITIES 

The Chaboche theory has been employed to predict a wide range of 
additional mechanical behavior of Inconel 718 at 1200°F, including: strain- 
rate and strain-rate-history effects, creep and relaxation behavior, load- 
unload- reload behavior under strain or stress control, and strain- 
controlled, fully- reversed cyclic behavior. 

Strain-Rate and Strain-Rate-Historv Effects 

Three different strain rates; namely, e ^ - 5x10 ^sec \ e ^ - 
5x10 ^sec \ - 5x10 ^sec \ have been arbitrarily considered in studying 
the strain- rate effect. The stress -strain curves at these rates are shown 
in Fig. 5, which shows a rate dependency of initial yield that could not be 
predicted by the original Chaboche theory. This figure shows that the 
extended Chaboche theory predicts a pronounced strain rate effect for the 
strain rates considered. It is also noted that the stress-strain curves 
tend to be closer to each other as e increases. This behavior is consistent 
with the experiments, which have shown that an asymptotic behavior can 
eventually be reached. 

A similar stress rate effect was predicted for three stress rates, 
obtained by multiplying the aforementioned strain rates by Young's modulus, 
as shown in Fig. 5. It is seen that the stress-controlled response 
overrides the strain-controlled response, the difference being attributed to 
the stress -strain nonlinearity of the constitutive equations. 

The response of viscoplastic materials is significantly influenced by 
their strain rate history [24,26]. Perhaps the most useful and widely used 
experiment developed to study this behavior is the incremental (jump) test, 
where a specimen is first prestrained at one strain rate and then the strain 
rate is abruptly changed to another value. A material is said to exhibit a 
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strain-rate-history effect if the jump (interrupted) and monotomic (pure) 
responses remain distinct after significant strains have been imposed [ 21 ]. 

Let the constant strain rates before and after a jump be denoted 
e ^ and respectively, and let 5 - € 2^ € l' r w here 6 > 1 indicates a strain 
rate increase, 0 < 5 < 1 corresponds to a strain rate decrease with 
increasing strain, and 5 < 0 implies strain rate reversal (unloading), Fig. 
6 . It can be shown that [27] 


* 

E 2 “ 


- ?< E - 


* 

V 


(38) 


where and are the slopes of the stress-strain curve immediately before 
and after the strain rate jump is imposed. If it is further assumed that E. 

« E, it follows that ~ E for 6 » 1. In other words, the initial 
response after the strain rate increment is imposed is elastic, as 
experiments show. In case of stress-rate jumps, the analogue of Eq. (38) is 


E 



(39) 


• • 

where 7 - a o/ a i t * ie ratio of the stress rates [27]. 

^ ^ *61 
Figure 7 shows the predicted monotonic responses at € - 10 sec and 

• -5 -1 1 

€. 2 - 5x10 sec as well as the corresponding interrupted responses for 

5-50 (low-to-high strain rate) and 5-0.02 (high-to-low strain rate) 
after a prestrain of 0.75%. 

For tensile loading, the saturation stress can be reached when each of 
the terms on the right hand side of Eq. (13) saturates. Numerical exercises 
with the Chaboche theory, however, have shown that R does not change 
substantially when small deformations, such as those produced during a 
tensile test, are considered, i.e., R » R q . Thus, the hardening process is 
essentially kinematic, and, in the presence of recovery, the saturation 
value of Y is smaller than a and can be determined by solving Eq . (32). 
This equation shows that the saturation value of Y is solely dependent on 
the strain rate and is, therefore, strain-rate-history independent. The 
difference observed in predicted monotonic and interrupted behavior is, thus 
caused by different values of R « R in the monotonic and jump tests. 
Since R q increases with e, the jump response lies below the uninterrupted 
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response for S > 1, i.e., the theory predicts a strain rate history effect 
typical of FCC materials [26]. Due to its treatment of R q as a constant, 
the original Chaboche theory cannot predict strain-rate-history effects. 

Prediction of the creep behavior of Inconel 718 at - 140 KSI, for 
load application rates of a ^ - 20 KSI /sec and a ^ ~ 0.20 KSI /sec, are 
presented in Fig. 8. In each case, the response from the moment of loading 
until the stress reached its hold value was predicted as well as the 
response for an additional 300 sec, with only the creep strains shown in 
Fig. 8. It was found that the creep strain corresponding to the higher 
stress rate overrides that corresponding to the lower rate. However, the 
total strain (not shown in Fig. 8) in the latter case was greater. This was 
expected, since at the lower stress rate, yielding initiated at a stress 
value much lower than - 140 KSI and by the time the stress value reached 
the hold value substantial hardening had already taken place. Thus, the 
material was capable of sustaining further deformations only at lower rates . 
These results suggest the importance of specifying the rate of load 
application, in addition to the hold stress, when reporting creep data. 

When steady-state creep is reached, Eq. (13) yields upon 
differentiation with respect to time 

Y - R - 0 , or Y — -R • (40) 

For tensile creep, Y is positive, and Eq. (40)^ is valid only if R is 
negative, i.e., when the material undergoes isotropic softening. Numerical 
solutions have shown that Inconel 718 at 1200°F starts to harden 
kinematically and soften isotropically with the initiation of creep 
deformations, with the rate of hardening being higher than the softening 
rate. However, a point is eventually reached where both rates are equal in 
magnitude. At this point Eq. (40)^ is satisfied and thereafter softening 
begins to dominate, thus allowing for increasing strain rates, as Fig. 9 
depicts for creep at 125 KSI. As mentioned earlier, there is experimental 
evidence that the increasing creep rates, at least in early stages, are not 
a manifestation of instability or impending failure [3,6]. Therefore, it is 
concluded that the Chaboche theory can predict realistic creep behavior, 
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including increasing creep rates following steady-state creep, caused by 
material softening. 

Relaxation Behavior 

In a stress relaxation situation, Eqs. ( 6 ) and (7) give 

a - - Ec” • (41) 

If the stress relaxes from a positive value, c M is initially positive and a 
is, therefore, negative. Experimental data on metals and alloys show that 
the stress normally relaxes at a decreasing rate until a steady state is 
eventually reached when a - 0 , i.e., when 

a - Y + R • (42) 

To examine the stress relaxation prediction of the Chaboche theory, a 
loading history was chosen, which consisted of straining at a rate of 5x10"^ 
sec to 1% strain. The strain was then held at this value for 50 minutes. 
The results are shown in fig. 10, and are consistent with the above 
discussion. 

Load-Unload-Reload Behavior 

A reversal in strain (stress) corresponds to an instantaneous jump in 
e(<j) with 6 ( 7 ) < 0. Three different cases can be distinguished (c.f. 

Fig. 6 ): (a) 6( 7 ) - -1, which corresponds to a reversal with constant strain 

(stress) rate magnitude, (b) 5 ( 7 ) « - 1 , which represents a large increase 
in the strain (stress) rate magnitude, and (c) 1 / 6 ( 7 ) « - 1 , which 
represents a large decrease in strain (stress) rate magnitude upon 
unloading. 

The loading-unloading behavior of Inconel 718 at 1200°F is now studied 
under strain (stress) control in a fashion similar to that of Ref. [28]. 
For strain control, a loading rate of - 5x10 ^ sec ^ was arbitrarily 

selected, along with unloading rates of - -5x10 ^ sec ^ (6^--10^«-l) , 

e 2 - -5x10 7 sec ' 1 (^ 2 — 1) . « 3 - -5xl0 ' 9 sec ' 1 (S^-10' 2 , or 1/Sj— 10 2 «-1) . 
In each case, unloading was initiated after a strain of 0.75% was reached, 
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where « 1, as was assumed in deriving Eq. (38). The results are shown 

in Fig. 11, where it is seen that the initial unloading behavior is 
inelastic and is sensitive to 5. For a large increase in the strain rate 

■jlf 

magnitude at reversal, the terms with 1/5 in Eq. (38) are negligible and E^ 

« E, i.e., the initial unloading response is elastic, which is in agreeement 

.2 

with the predicted behavior shown in Fig. 11 for 5 1 - -10 For 5* - -1, 

Eq. (38) gives E 2 - 2E, which means that the initial unloading response is 

inelastic. Figure 11 however, shows that the slope very rapidly 

approximated E as the stress decreased. For large decreases in strain 

rate magnitudes at reversal, the terms with 1/5 dominate and E 2 « E/|5|. 

Therefore, significant inelastic unloading response with a very large 

positive slope was predicted, resulting in near-vertical decay of the stress 

-2 

at the start of unloading. The inelastic strain increased when 5^ - -10 

was imposed, and can be thought of as stress relaxation, since the stress 

decreases at almost constant strain. 

For stress control, the stress rates were obtained by multiplying the 

prior strain rates by E. In each case, the loading period was adjusted so 

that a strain of 1% was obtained before unloading was started. The results 

* 

are shown in Fig. 12. For 7«-l, Eq. (39) yields a value of E 9 -E, i.e., the 

Z 2 

initial unloading response is elastic, as Fig. 12 depicts for 7 ^ - -10 . 

For 7 - - -1, the initial slope upon unloading can be found from Eq. (45) to 

be E 2 — -E^, i.e., the initial unloading response is inelastic with a 

negative slope that is equal in magnitude to the slope before unloading. As 

the stress decreases, however, the unloading response approaches a linear 

shape typical of elastic unloading, as shown in Fig. 12. For l/ 7«-l, Eq. 

(39) reduces to E 2 " - 1 nr | E^ . Consequently, the initial unloading slope is 

negative and much smaller in magnitude than the slope before unloading. The 

theory, therefore, predicts significant inelastic unloading behavior with 

small decreases in stress and very large increases in inelastic strain, as 

.2 

shown in Fig. 12, for 7 ^ - -10 . This latter case is analogous to a creep 
test, where the stress is held almost constant while the strain increases 
substantially. 

The analysis was further expanded by studying the effect of unloading 
rate on subsequent reloading behavior. Two cases were considered: (a) 
loading, unloading, and reloading at |f| - 5x10 ^sec , and (b) loading and 
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reloading at and unloading at - -5x10 ^ sec In each case, loading 

was applied until a strain of 0.75% was reached, unloading was initiated and 

maintained to a strain of 0.45%, and finally loading was continued to a 

strain of 1.0%. The results are shown in Fig. 13 for both cases. In the 

first case (6 - -1 at unloading), the initial slope upon unloading is 

positive and twice the elastic response, i.e., the initial response is 

inelastic. However, it rapidly approximates the linear elastic response 

and, therefore, the reloading response immediately approaches the monotonic 

_ 2 

response at the loading rate. In the second case (6 - -10 at unloading), 
however .significant inelastic deformations occurred upon unloading and the 
consequent reloading response at was considerably different from the 

monotonic response at the same rate. Note that differences in reloading 
behavior in the two cases are due essentially to the amount of inelastic 
deformations that occurred upon unloading. It is possible to conclude, 
therefore, that the Chaboche theory is capable of simulating a history- 
dependent reloading response that memorizes rate -dependent prior 
deformations . 

Cyclic Softening 

As mentioned earlier, Inconel 718 at 1200°F undergoes cyclic softening 

for approximately 10-20 % of its fatigue life. In Chaboche' s theory, 

cyclic softening is modeled by the isotropic hardening variable R, with its 

saturation value q smaller than its initial value R . The first few 

o 

hysteresis loops of a strain-controlled, fully- reversed cyclic test at |e| - 
4x10 ^ sec ^ and Ae - 2% were predicted and are shown in Fig. 14. In this 
figure, cyclic softening is manifested by a continuous decrease in the 
stress range and a corresponding increase in the inelastic strain range with 
cycling. 


V. CONCLUSIONS 

The modified Chaboche theory of viscoplasticity has been employed to 
model the viscoplastic behavior of Inconel 718 at 1200°F. A procedure has 
been developed and used for the determination of the material parameters of 
this theory from available uniaxial experimental data of standard mechanical 
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tests. The procedure for determining the parameters is not unique, but is 
considered appropriate for the available data, which were not adequate for 
the purpose of evaluating the Chaboche material parameters. 

The sensitivity study undertaken herein has resulted in an 'optimal' 
set of material parameters, which gave better overall agreement between 
theory and prediction. The study has also emphasized the influence of each 
material parameter on different aspects of material behavior. 

The predictive capabilities of the theory have been demonstrated for a 
variety of uniaxial loading conditions. The lack of adequate data base, 
however, precludes drawing detailed conclusions about the accuracy of the 
predicted behavior. Nonetheless, the Chaboche theory appears to offer a 
considerable promise for successfully modeling viscoplastic material 
behavior. 
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Table 1: Values of K , n , a, and c for Inconel 718 at 1200°F 


Ser , 
No, 

Ref. 

Test 

Strain 

Rate t 

Sec* 

No. of 
Data 
Points 

* 

K 

KSI 

★ 

n 

a 

KSI 

c 

1 

4 

cyclic 

8.333x10** 

4 

202.3 

0.1080 

24.46 

364.4 

2 

4 

cyclic 

3.333x10*, 

3 

188.7 

0.1100 

31.81 

347.0 

3 

4 

cyclic 

1. 667x10* f 

4 

156.6 

0.0920 

33.86 

349.0 

4 

5 

cyclic 

4. 000x10 *f 

12 

188.2 

0.1078 

31.49 

349.1 

5 

2 

cyclic 

4. 000x10* f, 
4. 000x10* ft 
2.667x10*;:' 

37 

172.2 

0 . 1017 

28.13 

355.0 

6 

2,5 

cyclic 

49 

187.9 

0.1096 

31.64 

347.5 

7 

9 

cyclic 

3 

214.9 

0.1562 

39.97 

301.5 

8 

4 

tensile 

1.670x10*^ 


166.9 

0.0540 

— 

— 

9 

4 

tensile 

3.330xl0* 5 


210.0 

0.0630 

— 

— 

^Some tests were performed at e 

- 10 ^ sec ^ and 

in*5 -1 

10 sec 




Table 

2: Calculation of R 

and p for the 

determination of b and q 

Cycle 

a 

€" t 

p 

Y 

R 

No. 

KSI 

% a 


KSI 

KSI 

i 

115.0 

0.530 

0.0265 

28.57 

39.04 

4 

109.8 

0.546 

0.0907 

28.71 

33.32 

8 

105.4 

0.561 

0.1568 

28.84 

28.83 

12 

102.4 

0.561 

0.2247 

28.95 

25.72 

16 

101.0 

0.588 

0.2942 

29.04 

24.25 

20 

99.42 

0.588 

0.3648 

29.04 

22.66 

24 

98.43 

0.588 

0.4354 

29.04 

21.67 

28 

97.43 

0.588 

0.5059 

29.04 

20.67 

32 

95.44 

0.591 

0.5766 

29.06 

18.66 

36 

95.44 

0.596 

0.6477 

29.09 

18.63 

40 

95.44 

0.596 

0.7192 

29.09 

18.63 


t 


Inelastic strain amplitude, Aa"/2 
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Table 3: Material parameters of Chaboche's theory, Inconel 718 at 1200°F 


Para- 

meter 

Description 

Initial 

Value 

Optimal 

Value 

E 

Young's modulus 

24.73x10 3 

24. 73X10' 3 

V 

Poisson’s ratio 

0.3356 

0.3356 

a 

Saturation value of kinematic 
hardening variable 

31.40 

30.00 

c 

Kinematic hardening exponent 

349 - 5 -10 
0.4x10 1U 

350 -° -10 
0.4x10 

7 

Coefficient of recovery 

m 

Recovery exponent 

6 . 942 | 

154.4 KSl n J sec 

7.000 r- 

155.0 KSI n Js< 

K 

Overstress parameter 

n 

Strain rate sensitivity parameter 

5.120 

5.100 

q 

Saturation value of isotropic 
hardening variable 

18.62 KSI 

50.00 KSI 

b 

Isotropic hardening exponent 

4.679 

3.75 

R 

o 

Initial, rate -dependent value of 
isotropic hardening 

Eq. (31) 

Eq . (37) 


Table 4: Effect of changing n by ±(3.5, 7.0)% on tensile behavior, 
Inconel 718 at 1200°F 


Strain 



Saturation Stress, KSI 


Rate. 

Exper . 



Predictions 


Sec 


0.930n 

0. 965n 

n 1.035n 

1 . 70n 

l.lxlO' 6 

134 

128.28 

129.20 

130.18 131.12 

132.33 

1.3x10* 5 

142 

143.57 

144.85 

146.23 147.62 

148.90 
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Table 5, Effect of changing R q by +(1.75, 3.5, 5.25, 7.0)% on creep behavior. Inconel 718 at 1200°F 
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Table 6: Effect of changing material parameters by values ranging from 




S~| W - 

w u cr 
< <o w 



63 




INELASTIC STRAIN 


Fig. 1. Illustration of two stress- 
inelastic strain curves at 
different strain rates. 



Fig. 2. Comparison of experimental 

and predicted mono tonic res- 
ponse at e = 1.333xl0 _5 sec _1 . 
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CREEP STRAIN 


0.450 




-1. 2-0.8 -0.4 0.0 0.4 0.8 1.2 

STRAIN, % 

Fig. 4. Predicted cyclic response at |e| = 
4xl0 -5 sec~ 1 and e = +1%; only ex- 
perimental peak stresses are avail- 
able. 
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STRAIN (STRESS) 


160 





00 4 

TIME 

CREEP TIME, sec 


Fig. 6. Illustration of strain (stress) 
rate changes. 

Fig. 8. Predicted creep response 
= 140 KSI for different 
rates. 

at 

loading 





CREEP STRAIN 



CREEP TIME, sec 


Fig. 9. Predicted creep response at - 125 KSI 
and o = 20 KSl/sec. 


Predicted load-unload res- 
ponse at different unload- 
ing strain rates. 



0.0 800 1,600 2,400 3,200 0.0 0.8 1.6 2.4 


TIME, sec 


STRAIN, X 


Fig. 10. Predicted relaxation response at 
e = 1% 


Fig. 12. Predicted load-unload res- 
ponse at different unload- 
ing stress rates. 
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STRAIN, X 

Fig, 13. Predicted load-unload- 
reload response at dif- 
ferent unloading strain 
rates . 



STRAIN, X 

Fig. 14. Predicted cyclic response 
at | e | = 4x10 5 sec 1 and 
e = + 1% . 
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